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Abstract 

We study the elastic response of bilayer membranes with fixed projected area to both stretching 
and shape deformations. A surface tension is associated to each of these deformations. By using 
model amphiphilic membranes and computer simulations, we are able to observe both the types of 
deformation, and thus, both the surface tensions, related to each type of deformation, are measured 
for the same system. These surface tensions are found to assume different values in the same bilayer 
membrane: in particular they vanish for different values of the projected area. We introduce a 
simple theory which relates the two quantities and successfully apply it to the data obtained with 
computer simulations. 
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Bilayer membranes are composed of amphiphilic molecules whose hydrophobic part is 
strongly insoluble in water. Such molecules tends to form interfaces with the solvent in 
order to reduce the interfacial energy a bilayer is thus formed by two adjacent sheets 
of amphiphilic molecules separating two aqueous phases. The interest in such structures 
resides in the fact that they play an outstanding role in the organization of biological cells: 
bilayer membranes of lipids form the basic structure of a cell's plasma membrane and of 
internal membranes, which surround the organelles in eukaryotic cells, such as the nucleus, 
mitochondria and the Golgi apparatus 2], |J. 

In a fluid membrane fluctuating freely in a solvent, all internal degrees of freedom, related, 
e.g, to the hydrocarbon chain conformation, and to to the local molecular density, relax on 
a fast time scale. The membrane is then characterized by a zero shear modulus, i.e. any 
shear deformation induces a flow of the amphiphilic molecules within the membrane. Such a 
system is then subject to only two types of elastic deformations: bending deformations and 
stretching deformations. The first are deformations normal to the membrane plane which 
change the membrane shapes, while the latter are in-plane deformations which change the 
local area per molecule projected onto the membrane midplane. 

The surface of a fluctuating membrane can thus be viewed as an interface with the solvent, 
and we will call fluctuation surface tension the surface tension of the bilayer membrane 
associated to shape fluctuations. The surface tension of an interface saturated by surfactant 
molecules is expected to vanish (see, e.g., Q), and thus the fluctuation surface tension of 
a fluid membrane fluctuating freely in a solvent is usually assumed to be zero. However in 
the case of a constrained bilayer membrane, such as those fixed on a frame, the geometrical 
constrain can lead to a non vanishing surface tension. Non vanishing fluctuation surface 
tension have been observed, e.g., in [5|. 

On a macroscopic length scale, a membrane fixed on a frame can be treated as an elastic 
sheet: if one of the sides of the frame is movable and one tries to change the frame area by 
pulling or pushing on this side, an elastic force will appear which will tend to restore the 
preferred value of the frame area. The parameter describing such elastic behaviour is the 
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mechanical surface tension and can assume non negative values too 

Intuitively, one would expect that the two surface tensions associated with the two differ- 
ent membrane deformations, namely shape deformation and lateral stretching, are different 

n 

quantities, see ref. |£| and discussion therein. However, in experiments, they are usually 
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measured by observing the one or the other type of deformation, and so one cannot compare 
them experimentally for the same system. On the contrary, the two surface tensions can 
be measured using computer simulations of model membranes, and this has been done in 
two prev.ous wo rk s Q fi. But in one case, the two surface tensions have been found to 
be equal within the statistical errors P|, while in the other case they have been found to be 
proportional with respect to each other [lo| . 

The aim of this paper is thus the systematical study and comparison between the fluctu- 
ation surface tension and the mechanical surface tension, measured by observing both the 
types of deformation in the same bilayer membrane. In the present paper we consider bilayer 
membrane with fixed projected area. The projected area correspond to the area of the frame 
in the picture drawn above. By using model bilayer membrane and computer simulations, 
we measure both the surface tension associated to shape deformation and that associated 
to stretching deformations. We show, for the first time, that in the same system these two 
quantities assume different values, and in particular they do not vanish for the same value 
of the projected area. The model used here for the computer simulations, have been largely 
used in previous works and has been quite successful in describing many dynamic and elastic 
properties of real bilayer membranes |_6 
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The paper is organized as follows: in section |1J we describe the model bilayer, together 
with the computer simulation techniques which we use in the present work. In section ITT1 we 
discuss the theory of the mechanical surface tension and the method adopted to measure 
it. We then show the results for such quantity obtained by computer simulations. In 
section UTT] after reviewing the classical elasticity theory for the shape fluctuations of bilayer 
membranes, we show the results for the shape fluctuation surface tension, obtained by 
computer simulations. In section IIV1 the two quantities are compared, and a simple theory 
which connects them is introduced and discussed. We discuss our results and conclude in 
section E| 

In the following, we consider a bilayer membrane made up of N amphiphilic molecules, 
at temperature T, spanning a square frame of area A p = L 2 , its effective area being A. 
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I. MODEL BILAYERS AND COMPUTER SIMULATIONS 



A. Coarse grained model of amphiphilic bilayers 
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We adopt the same coarse grained model introduced and used in refs. 
turned out to be an effective model to study several properties of bilayer membranes such 
as surface tension, bending rigidity and diffusion characteristics. The amphiphilic molecules 
are represented as linear chains of beads, a single bead representing the molecule hydrophilic 
head (H), or several CH2/CH3 groups of the amphiphile hydrocarbon chain (C) see fig. 
The water molecule is also represented as a single bead (W). The hydrophobic interaction 
of C with W and H particles is modelled by soft core potential 
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Uscinj) = 4e 
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(1) 



while attractive interactions W-W, W-H, H-H, C-C, are modelled by a Lennard- Jones po- 
tential 

W^4 £wJ [(i-) I2 -(i-) 6 ], (2) 

where a//3e {W, H, C}. Adjacent beads along a single model molecule interact via the 
harmonic potential 

= k 2 ( r v+i - , (3) 

where i and i + 1 indicate two successive particles along the chain. The effects of hydrocarbon 
chain stiffness is modelled by letting all the particles within a single amphiphile molecule 
interact via the three-body bending potential 



^3( r i-l,i> r i,i+l) — ^3 

where r i i+1 = r i+1 - r», see fig. EJ 
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(4) 



B. Simulation method and model parameters 

In the present work we combine both Monte-Carlo (MC) and Molecular Dynamics (MD) 
simulation methods. Our simulations are carried out for a cuboidal boxes of constant volume 
and periodic boundary conditions. The MC algorithm is used to make the system relax 
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towards configurations of minimal potential energy, the output of 500 x N MC steps are 
used as starting configurations for the MD part of the simulations. 

For the parameters characterizing the above potentials we choose the values used in 
refs. B y Q Q. F„ r the faction ranges of the U potential we take , = 1/3 n.n 
which is of the same order as the LJ lengths of interactions for pairs of CH 2 groups, CH 3 
groups or water molecules, as discussed in The characteristic length of the repulsive 
SC potential (jlj), is taken to be £' = 1.051 as in with this choice the hard-core repulsion 
of the SC potential is of the same order of the repulsive part of the LJ potential. We take 
e = 2/N Av kj (AV is the Avogadro number): this value is bigger than the LJ energ.es 
for pairs of CH 2 and/or water molecules reported in [14j . but it takes into account that in 
our model one C particle corresponds to three or four CH 2 groups The strength of the 
harmonic bond potential @3) and of the three-body bending potentials are taken to be 
k 2 = 5000e/£ 2 and k 3 = 2e p. 

Since we run MD simulations, two additional parameters are involved in the present 
model: the masses of the different beads, which enter the equations of motions, and the 
time step 5t, used to discretize such equations. For simplicity's sake, all the beads are taken 
to have the same mass m = 0.036/ N^ v kg, which is approximately the mean value of the 
mass of a water molecule and the mass of four CH 2 groups. Using this value of m one 



obtains the characteristic time scale r = y^m£ 2 /e =1.4 ps. In the MD simulations, the 
integration time step St is taken to be 5t = r/2000 = 0.7 fs. The leap-frog algorithm at 
constant temperature |l^ with /c^T = 1.35e (corresponding to T = 325 K) is used for the 
integration of the equations of motions. In the following, if not differently specified, all the 
quantities will be expressed in units of the tree basic parameters £, e and m. 

In the present paper we consider bilayers with three different values of the number of 
amphiphilic molecules iV = 512, iV = 768 and iV = 1152. For each value of N, the 
simulations are run with different values of the projected area A p , which corresponds to 
the area of the simulation box side parallel to the bilayer. In changing the projected area, 
we keep constant the number of water particles Nyy, as well as the overall simulation box 
volume V, i.e., we keep the overall system density constant. The values of N, Ny/ and V 
used in the simulations described in the following are reported in table |U 
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N 


N w 


V(£ 3 ) 


512 


3200 


8640 


768 


4800 


12960 


1152 


7200 


19440 



Table I: Number of amphiphiles N, number of water particles Nyy and simulation box volume V 
in I s units, used in the simulations described in the text. 

II. STRESS TENSOR AND MECHANICAL SURFACE TENSION 



On a macroscopic length scale, a membrane fixed on a frame can be viewed as an elastic 
sheet: the system will be characterized by some equilibrium value of the frame area where 
the force exerted on the frame sizes vanishes, and small changes in the frame area cause 
restoring forces to appear. If F is the free energy of a framed bilayer, the surface tension E 
conjugated to its projected area A p is defined as 

„ dF 



dA r . 



(5) 



In the present work, the bilayer surface tension is measured using a method first intro- 
duced by Schofield and Henderson [l^, and then extended by Goetz and Lipowsky A 
fluctuating membranes can be considered isotropic and homogeneous along any direction 
parallel to the membrane plane, if the amplitude of the fluctuations is not too high. Thus, 
the system stress tensor will be a function of the coordinate z perpendicular to the bilayer 
plane. The surface tension is related to the system stress tensor via 



(6) 



where ^t(z) and ^n(z) are the components of the stress tensor perpendicular and normal to 
the bilayer surface, respectively. The integral can be extended to infinity because the stress 
tensor in isotropic in the bulk water and so T>t(z) = Hn(z). By noting that the stress tensor 
is defined as the negative of the pressure tensor, one sees that eq. (JBJ) is equivalent to the 
original expression of the surface tension of a planar liquid- vapor interface, as formulated 
by Kirkwood and Buff [16]. The macroscopic stress tensor I] can be expressed in terms of 
the microscopic stress s which depends on the momenta and on the positions of the system 
particles within a small volume. The microscopic and the macroscopic stress tensor are 



related via 

S = (s) , (7) 

where the brackets denote ensemble average. The microscopic stress tensor s consists of two 
contributions one arising from the kinetic energy and the other from the interaction energy 
of the system particles within a small volume s = s^- + s in . The kinetic contribution to 
the macroscopic stress tensor, as given by eq. (jZJ), can be neglected, since it is isotropic, and 
must vanish on average. Thus, one is left with the interaction part of the microscopic stress 
tensor s in only, and eq. (0) becomes 

£ = <s> m , (8) 

Furthermore a planar bilayer membrane in solvent can be considered homogeneous along any 
direction parallel to the the bilayer plane, and so is the stress tensor. The system can thus 
be divided in thin slices parallel to the bilayer plane, and the stress tensor can be averaged 
over each of these slices. Thus the microscopic stress tensor is given by Jf| Q, 

s in (z) = f(zi, Zj, z)-£y Yl ® r tf > ( 9 ) 

i>j 

where the sum is extended to all the system particle pairs, AV is the volume of the slice, the 
function f(zi,Zj,z) determines the actual contribution of the pair i,j to the stress tensor 
in the slice of coordinate z, and the symbol £g> denotes the tensorial product. If both the 
particles are inside the current slice we take f(zi,Zj,z) = 1 Let Az be the thickness 
of each slice, if both the particles are external to the current slice and on opposite sides, 
we take f(zi,Zj,z) = Az/\zi — Zj\, where, using periodic boundary conditions, the shortest 
distance \z{ — Zj\ between them is considered. If just one of the two particles lies within the 
current slice, we take f(zi,Zj,z) = Az/(2\zi — Zj\). In all the other cases, the contribution 
of the particles i,j to the stress tensor associated to the slice of coordinate z vanishes, and 
thus we take f(zi,Zj,z) = 0. Inspection of eq. suggests that its rhs can be viewed as 
a local density of the macroscopic virial tensor W = Yli>j ^ij ® r i?> an d thus the function 
f(zi,Zj,z) determines the fraction of virial tensor density to be associated to the slice of 
coordinate z. 



t In ref.jg, the expression for the contribution to the stress tensor of the three body interaction (@J was 
sligthly different. However, a three body interaction can be expressed as the sum of pairwise interactions 
|l5| , and thus eq. 10 is correct also for the three body potential . 



A. Measurements of the microscopic stress tensor 



Using the model amphiphilic bilayer described in section U the microscopic stress tensor 
Si n (z), as given by eq. (|9jL was measured for the three values of the number of amphiphilic 
molecules N here considered, N = 512, N = 768, and iV = 1152. For each value of N, s in (z) 
was measured for different values of the bilayer projected area A p keeping constant the total 
volume V, see table [I] The simulation boxes have been divided into 140 ^-slices, which 
corresponds to Az ~ 0.1£, for the simulation box dimensions here used. The off-diagonal 
elements of the system stress tensor are found to vanish as expected, as well as the average 
of its kinetic part (data not shown). For each value of A p , the mechanical surface tension E 
is obtained using eq. (jUJ), by performing a discrete integration. For each value of N and A p 
we run 5 simulations of 2 x 10 5 MD steps, and a run of 6 x 10 6 MC steps was interposed 
after each run of 2 x 10 5 MD steps. The MC steps are inserted to allow the system to reach 
regions of its phase space not easily accessible with a single MD trajectory. 

It is worth noting that the code for the stress tensor measuring is extremely time- 
consuming: a single run of 2 x 10 5 MD steps, for a system with N = 1152 amphiphiles, needs 
about one month of machine time, on a computer equipped with a 1 GHz processor. The 
relative fluctuations of the surface tension with respect to its average value are very large, 
as found in other work |(|. Following ref. Q], in order to reduce the effect of short time 
fluctuations, the total simulation time of each run is divided in blocks of 5000 MD steps, 
and the surface tension is subaveraged over each of these blocks. This average value over 
the 5000 MD steps is then used as the sample value for the current block. 

The surface tension obtained from simulations is plotted in fig. |21 as a function of the 
projected area per amphiphile a p = A p /(N/2), for the three values of N here considered. 

According to the classical elasticity theory, the parameter describing the effective com- 
pressibility of an elastic sheet to a change in its projected area is the area compressibility, 
Ka defined as 



K A = A p 



<9£ 



(10) 



9A„ 



By integrating this last equation, one obtains 




(11) 



where the last equality holds for values of A p close to A p , which is the area at which £ 
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vanishes. Upon integration of eq. (jllj) . we obtain the classical expression for the free energy 
of a stretched (or compressed) elastic sheet as a a function of A p around the equilibrium 
value A p 

*- = g(A>-4) 2 . (12) 

where the reference free energy at A = A p has been chosen to be equal to zero. Inspection of 
figure El indicates that the surface tension S obtained by simulations follows the behaviour 
predicted by eq. (jllj) : this quantity is linear within a range of values around the equilibrium 
area A p . We find that for A p <C AL (a p ~ 2.065) the surface tension £ is no longer a 
monotonous function of a p : this is probably due to the fact that, for such small projected 
area per molecule, the system exhibits buckling, and the membrane cannot be considered 
flat on average. In this case, the basic hypothesis that the membrane is isotropic and 
homogeneous along the the z-axis is not fulfilled. However, such hypothesis is required for 
measuring £ with the method discussed in this section, and thus the results obtained for 
very small values of a p must be inaccurate. Thus in the following we will consider values of 
a p such that S is a monotonic function of a p , i.e. a p > 2.065. The values of the equilibrium 
area, of the area compressibility and of the slope of £ as a function of a p , for the three system 
sizes here considered, are listed in table |H] In the same table, we also report the results for 
a smaller system which was considered in a previous work [6] . The results listed in table |H] 



N 


4 (^ 2 ) 


K A (e/t 2 ) 


dZ/da p (e/£ 3 ) 


128 


2.15 ±0.02 


11.8 ± 1.5 


5.5 ±0.7 


512 


2.12 ±0.01 


11.7 ± 1.5 


5.5 ±0.7 


768 


2.12 ±0.01 


9.4 ±1.1 


4.5 ±0.5 


1152 


2.12 ±0.01 


8.8 ±0.8 


4.1 ±0.4 



Table II: Equilibrium projected area per amphiphile compressibility Ka for different 

system sizes, and slope of £ with respect to the projected area per amphiphile a p . The values for 
N = 512,768, 1152 are obtained by linear fit of the data shown in fig. [2] to ea. l|llj) . in the range 
where X is a linear function of a p . The data for N = 128 are taken by ref. j^|. 

suggests that the equilibrium area A| does not depend on the system size, but is an intrinsic 
properties of the model amphiphilic molecule used here. However, the area compressibility 

n 

decreases as a function of the system size, as already found in p : this decrease is due to the 
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fact that a larger number of oscillation modes are introduced in the system as its projected 
area is increased. Thus, the bilayer becomes easier to compress or to stretch if its total 
projected area A p increases while the projected area per amphiphile a p is kept constant. 
Note that the values of the area compressibility found for the present model range between 
263 and 360 mJ/m 2 , see table |HI Such values are slightly greater than those observed for 
real lipid bilayers, that are in the range 140-240 mJ/m 2 [l8j, but are very similar to those 
values found m other works on computer simulations of amphiphilic bilayer with atomic 
resolution k3, 



III. ELASTICITY HAMILTONIAN AND THE FLUCTUATION SURFACE TEN- 
SION. 

In the Monge representation, the shape of the bilayer is described by the height function 



h(x, y) which measures the distance of its midsurface from the reference (x, y) plane, 
classical elasticity Hamiltonian of fluctuating membranes, for small fluctuations, reads 



\a{Vhf + l -K{V 2 h) 2 



^he 





H = aA p + / dxdy 

J A p 

= °A> + ^2 J2 W + K ? 4 1 l^l 2 ' ( 13 ) 

q 



where /i q define the Fourier coefficients of h(x, y) 

frq = jr J dxdy exp [ir ■ q] h(x, y). (14) 

The two parameters appearing in eq. (|T3*j) are the bending rigidity k, describing the resistance 
of the system to bending, and the surface tension a, which takes into account the contribution 
of the bilayer total area A to the system energy. In fact, eq. (|13J) can be rewritten as 

H = aA + j ' dxdy^K(V 2 h) 2 . (15) 

The membrane fluctuation spectrum, defined as S(q) = \h n \ 2 , depends only on q = |q| and 
exhibits the functional form (2]] 

S ^ = A i 4V ( 16 ) 

Ap {crq 2 + fcq ) 
10 



Equation (|16p has been obtained under the hypothesis that the bilayer can be modelled as 
a geometrical surface, neglecting its thickness. Thus eq. (j!6|) holds only for wavenumbers 
smaller than an upper bond q*, which is usually taken to be q* ~ 2n/do, where do is the 
bilayer thickness. We find that the model bilayer thickness is do ~ 6£, independently of the 
system size. For q > q*, the shape of the bilayer is characterized by the local protrusions of 
the amphiphiles, and thus the fluctuation spectrum is dominated by a local surface tension 
term, see refs 0[llI[l3|. However these sub-optical modes are not accessible in experiments, 
and thus, in the present work, we will focus on the small wavenumber, and on the effective 
values of n and cr that can be obtained by eq. (fT7)|) . By sampling the height field h(x,y) 
during the simulations, one can obtain the mean fluctuation spectrum (fTtjj) . By fitting S(q) 
to the rhs of eq. (JTHJ), for q < q*, the values of the bending rigidity k and of the surface 
tension a can thus be estimated. The fluctuation spectrum shape depends on the projected 
area: and thus both k and a depend on this parameter. By changing the value of A p one 
can achieve the tensionless (a = 0) state Q] . The surface tension has been measured for 
the three system sizes here considered iV = 521, 768, 1152, for different projected area, and 
keeping constant the overall volume, see table HI For each value of the projected area, we 
run 10 x 10 5 MD steps: a run of 6 x 10 6 MC steps has been interposed after each run of 10 5 
steps. The fluctuation spectrum has been sampled during the MD simulations, every 5000 
MD steps. In figJ3]the fluctuation spectrum of the larger system here considered is plotted 
as a function of the wavenumber q for different value of the projected area per amphiphile 
a p : the surface tension a vanishes for a given values of a p = a p , which depends on the 
microscopic details of the system (on the model parameters in the case of simulations), 
while for a p > a* (a p < a*) it is greater (smaller) than zero. This is clearly shown in figure 
where the surface tension a is plotted as a function of a p for the three values of N here 
considered. 

By linear fit of the data in figure 01 we obtain the tensionless area a* and the slope of 
the surface tension a as a function of a p , the results are listed in table ITTT1 As in the case of 
the mechanical surface tension E, the slope of the curve a(a p ) decreases as N is increased. 
Comparison of tables |H] and IIHI indicates that the slopes of E and a as functions of a p are 



different for a given system size. 
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N 


4 (? 2 ) 


da/da p (e/£ 3 ) 


512 


2.085 ± 0.005 


5.4 ±0.7 


768 


2.09 ±0.01 


4.0 ±0.4 


1152 


2.09 ±0.01 


3.6 ±0.4 



Table III: Equilibrium projected area per amphiphile a*, and slope of a with respect to the projected 
area per amphiphile a p . The values for the three system size are obtained by linear fit of the data 
shown in fig. EI 

IV. COMPARISON BETWEEN THE TWO SURFACE TENSION a AND £ 

Comparison of figure|3]with figure 121 shows clearly that the projected area per amphiphile 
a* ~ 2.09 where the surface tension a vanishes is different from the projected area per 
amphiphile a p ~ 2.12, at which the mechanical tension £ vanishes. Since these two quantities 
have been measured independently, for each value of the projected area, this last result 
suggests that one is dealing with two different quantities. In addition, the two vanishing 
areas a* and ap do not exhibit any dependence on the system size N, and thus one can argue 
that they are two independent intrinsic properties of the amphiphilic molecule. 

Comparison of tables ITT1 and ITTT1 and inspection of figure clearly indicate that the curves 
S(a p ) and a(a p ) exhibit different slopes. Furthermore, if one plots on the same figure a and 
£ as functions of a p , in the range where £ is linear with respect to a p , the two data sets 
appear clearly shifted the one respect to the other, see figure where these two quantities 
have been plotted for the N = 1152 case. 

In order to relate the two surface tension a and E, we need first to define the statistical 
ensemble which characterizes the system considered here, i.e., the minimal set of thermo- 
dynamical variables which fully describe the system state. The classical argument that the 
system area readjusts to its preferred value (see, e.g., jj|) cannot be invoked here because 
of the constrain represented by the frame area (simulation box area). The effective bilayer 
area A can be evaluated using a discrete version of the formula 

A= f dxdyy / l±(V/i) 2 , (17) 

J A p 

where the field h(x, y) is sampled during the MD simulations. The effective area per am- 
phiphile a is plotted in fig. |H1 as a function of a p . Inspection of this figure clearly indicates 
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that the effective area a is not independent of the projected area a p , but is rather a function 
of a p . Thus, the thermodynamical ensemble characterizing the framed bilayer here consid- 
ered is the (T, A p = L 2 , N) ensemble, where L is the size of the simulation box face parallel 
to the bilayer. 

The free energy in the (T, A p = L 2 , N) ensemble for a bilayer membrane system is 



F = -k B T\nZ = -k B T\n (J Vhe- pHm) \ 

/N-l 

~ -k B T In I 11 



i=0 



dhj 



-k B T\n\T[ J 



A 

dh 



<! e -fSH({h q }) 



aA p — k B T In 



crA p — k B T In 



n 

q 

n 



aA p + \rk B T In 



Na p \ 



dh. 



Na p \ 



2tt 



X 2 a p /3(aq 2 + nq 4 ) 

X 2 a p (aq 2 + nq 4 ) 
2nk B T 



(18) 



The parameter A has the dimension of a length and is the analogous of the De Broglie wave- 
length in the ideal gas partition function. Without introducing this parameter, the partition 
function would be dimensionful. The equation of state for the intrinsic area is 



(A) 



WI4> =Ar+ h BTE _±_ . 
da 2 ^ a + nq z 



(19) 



which is a well known result, see, e.g, [22j. Let us now introduce the quantity N', which 
represents the number of membrane patches which fluctuate independently (N f < N/2): 
the summations in equations ()18j) and (|19jl run over these N' wavemodes, and thus we have 
= N'. Note that, in the large A p = L 2 limit, N' can be estimated as follows 



N' 



( L Y 



q 

L 2 



T 2 

dq 2 = — 
A, 



qdq 



-^{q -q mi n)^—, 



(20) 



where we have taken q* ~ 1 (see discussion is section ITTT|). and q* ^> q m i n — 2tt/L ~ 0, in 



the large system size limit. Taking into account that q = (2ir) (n 2 , + n 2 )/L , and since the 
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surface tension E is the thermodynamic conjugate of A p = L 2 , as defined by eq. ©, we have 

dF 1 k B T ^ nq* 

~dA p ~ ° ~ 2^47^ Kq 4 + aq 2 



AT,A P ) = — = g --_ > (21) 



q 



2 A p ^ nq 2 + a 2 A p 1 1 

= W a _kBTV j (23) 
2 Ap 

where we substituted eq. ()19|) into eq. (j22J). Equation relates thus the frame surface 
tension E with two of the variables which define the actual ensemble T, A p , and with the 



fluctuation surface tension o. Comparison of tables [TT1 and IIHI indicates that, for a given 
system size, the slope of E(a p ) is larger than the slope of cr(a p ): since A > A p by definition 
of intrinsic area, the first term on the rhs of eq. (|23j) accounts for the difference of the 
slope between the two curves. Note that the second term on the rhs of eq. (|23j) takes into 
account the effect of the entropic elasticity on the system surface tension. This term is 
proportional to the system temperature, and inversely proportional to the system effective 
area per amphiphile a' p = A p /N': the analogy with the case of polymers can be immediately 
drawn. It is worth noting that a similar results was obtained by Farago and Pincus in 
ref. (23}], where a continuous expression for the free energy (JTSJl was used, and where the 
surface tension a was taken to depend explicitly on the effective area A. A first rapid check 
of equation can be done by noting that it predicts that, if the surface tension a vanishes, 
the frame surface tension E is negative. This is confirmed by inspection of figures EJ and 0] 

In equation ()23|) the only adjustable parameter is N', while a can be estimated as de- 
scribed in section IIHI and A can be sampled during the MD simulations using eq. (JT7j) . 
One can thus compare the values for E predicted by eq.(J23J) with those directly measured 
as described in section |H] In the following, the optimal value of the parameter N' will be 
determined by fitting eq. (j2S| to the measured values plotted in fig. El for the three values 
of N here considered. With this fit we find N' = 155 for N = 512, N' = 181 for N = 768 
and N' = 213 for iV = 1152. Using this values for N', we plot the measured and the cal- 
culated values for E, as a function of a p , for iV = 512 fig. [?l for iV = 768 fig. |HJ and for 
iV = 1152 fig. El Inspection of these figures indicates a good agreement between the values 
of E predicted by eq. (j2*3|) and those calculated as described in section |TT1 

We now consider the scaling behaviour of the quantity N' as a function of the system 
size. In fig. [TUl N' is plotted as a function of the tensionless projected area A* (projected 
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area at which a = 0), and as a function of the number of molecules N. Inspection of this 
figure, suggests that N' is a linear function of A*, and the slope of such function is in good 
agreement with that predicted by eq. (|20|). It is worth noting that, as discussed is section 
11111 the projected area per molecule a* = A*/(N/2) is independent of the system size, and 
thus we have N' ~ A* J (Att) ~ N, as one would expect in the large N (or large A p ) limit. 

Note that the intercept of the line plotted in fig. El is non-zero, while eq. ()20|) predicts a 
vanishing intercept in the large projected area limit. However the argument used to obtain 
eq. ()2U|) is no longer valid for small value of A p (N), and thus in this limit such an equation 
is incorrect. 



V. CONCLUSION 



In the present paper we have measured both the surface tension a appearing in the 
Hamiltonian which governs the bilayer shape fluctuations Qlfljl . and the surface tension E 
which governs the elastic response of the system to a change in its projected area A v . Using 
computer simulations, we have measured independently these two quantities for bilayer with 
different projected area. Our results indicate that the two surface tensions have different 
values for a given value of A p . In particular the two projected area per amphiphile a* and a P 
are different, where a* is the projected area per amphiphile at which the surface tension a 
vanishes, while a| is the projected area per amphiphile where the surface tension £ vanishes. 
Furthermore, the two surface tensions are found to exhibit different slopes as functions of 
the projected area per molecule a p , for a given system size. Using a simple thermodynamic 
argument, we manage to relate the two quantities a and E, and the relation between them 
which we found, eq. ()23|) . nicely fits the data that we obtain from simulations. This is the 
most important result of the present work: it indicates that eq. ()23|) succeeds to capture the 
basic relation between the two quantities. Note that in a previous work Q| , where computer 
simulations of model membrane were considered, the two surface tensions were found to be 
proportional the one respect to the other, but no argument was introduced to explain such 
an effect. 

The difference between the two quantities, does not appear to be due to some size effect, 
since the vanishing areas a* and dp do not depend on the system size. The scaling behaviour 
of the effective wavemode number N' as a function of the system size is a rather delicate 
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issue. In the present work we find that N' scales linearly with the number of molecule N 
(the tensionless area A*), but in a restricted range of values of N (A p ). We were limited in 
the choice of the values of iV by the considerable amount of computation time required by 
the simulations. For larger values of N one would expect that N' grows either as iV or more 
slowly than N. If this were the case, since A p ~ N, the second term of the rhs of eq. (|23|) 
would vanish in the large N limit. In this case eq. (|23p would read £ ~ A/A p ■ a. We plan 
to consider larger system size in a future work. 

In conclusion, the results contained in this paper puts in evidence that in a bilayer mem- 
brane, the surface tension a appearing in the elasticity Hamiltonian (jlHjl and the mechanical 
surface tension £ are two distinct thermodynamic quantities, which have to be measured 
independently for a full characterization of the system elastic properties, although for a 
constrained system, the geometrical constraints impose a relation between them. 

It would be interesting to measure these two quantities in a real system, e.g, a lipid vesicle 
manipulated with a micropipette. On the basis of the results presented in this paper, one 
would expect that the geometrical constrains imposed by the micropipette would lead to a 
difference between the two surface tensions. 

Acknowledgments 

The author is indebted to R. Lipowsky for introducing him to this topic, for many in- 
teresting discussions and for a critical reading of this manuscript. I am grateful to J.B. 
Fournier for the long and interesting discussions on the two tension issue. I also thank L. 
Peliti for many discussions, and for his encouragements. This work was partially supported 
by MIUR-PRIN 2004. Finally, I thank the CRdC AMRA for the use of its computational 
resources. 



[1] J. Israelachvili, Intermolecular and Surface Forces, (Academic Press, London, 1992). 

[2] B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts, and P. Walter, Molecular Biology of 

the Cell, (Garland, New York, 2002). 
[3] R. Lipowsky, Vesicles and biomembranes in Encyclopedia of Applied Physics Vol. 23, pages 

199-222, (VCH Publishers, Weinheim and New York 1998). 

16 



[4] Safran S. A., Statistical Thermodynamics of Surfaces, Interfaces and Membranes (Addison- 

Wesley, Reading, Mass. 1994). 
[5] S.J. Marrink and A.E. Mark, J. Phys. Chem. B 105:6122-6127 (2001). 
[6] R. Goetz and R. Lipowsky, J. Chem. Phys. 108:7397-7409 (1998). 
[7] E. Lindhal and O. Edholm, J. Chem. Phys. 113:3882 (20 00). 
[8] F. David and S. Leibler, J. Phys. France 1:959-976 (1991). 
[9] O. Farago and P. Pincus, J. Chem. Phys. 120: 2934-2950 (2004). 
[10] J. Stecki, J. Chem. Phys. 120: 3508-3516 (2004). 

[11] R. Goetz, G. Gompper and R. Lipowsky Phys. Rev. Lett. 82:221-224 (1999). 

[12] A. Imparato, J.C. Shillcock and R. Lipowsky Eur. Phys. J. E 11:21-28 (2003). 

[13] A. Imparato, J.C. Shillcock and R. Lipowsky Europhys. Lett. , 69: 650-656 (2005). 

[14] E. Egberts and H.J.C. Berendsen J. Chem. Phys. 89:3718-3732 (1988). 

[15] M.P. Allen and D.J. Tildesley, Computer Simulations of Liquids (Oxford, Oxford, 1987). 

[16] J.G. Kirkwood and F.P. Buff, J. Chem. Phys., 17:338 (1 949). 

[17] P. Schofield , J.R. Henderson, Proc. R. Soc Lond. A MAT 379 ( 1776): 231-246 1982. 

[18] D. Marsh , CRC handbook of lipid bilayers, Boca Raton, CRC Press, (1990). 

[19] P.B. Canham, J. Theor. Biol 26:61-81 (1970). 

[20] W. Helfrich, Z. Naturforsch., 28c:693 (1973). 

[21] J.F. Lennon and F. Brochrard, J. Phys. (Paris) 36, 1035 (1975). 

[22] D. Roux, Physica A 172, 242-257 (1991). 

[23] O.Farago and P. Pincus, Eur. Phys. J. E, 11:399 (2003). 



17 



t + 1 



Figure 1: Cartoon of the model amphiphilc molecule used in this paper. 
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Figure 2: Plot of E as a function of the projected area per amphiphile a p for the three values of N 
here considered, N = 512, N = 768 and iV = 1152. 
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Figure 3: Fluctuation spectrum 5 as a function of the wavenumber q, for different values of the 
projected area per amphiphile a p , for bilayers with N = 1152. The tensionless state (a = 0) 
corresponds to the projected area per amphiphile a* = 2.095. The dashed line is obtained by 
fitting the fluctuation spectrum of the tensionless system to eq. ()16|) . with a = 0, for q < 1. 
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Figure 4: Plot of a as a function of the projected area per amphiphile a p for the three values of N 
here considered. 
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Figure 5: Plot of the fluctuation surface tension a and of the mechanical surface tension £ as 
functions of the projected area per amphiphile a p , for the N = 1152 system. The lines are linear 
fits to the two sets of data. The two sets of data appear to be clearly shifted and tilted with respect 
to each other. 
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Figure 6: Plot of the measured effective area per amphiphile a, as a function of the projected area 
per amphiphile a p (simulation box area), for the three values of N here considered. 
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Figure 7: Plot of the measured surface tension £ and of the estimate of the same quantity, as given 
by eq. (|23|). as functions of a p , for a system with N = 512 amphiphiles. 
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Figure 8: Plot of the measured surface tension £ and of the estimate of the same quantity, as given 
by eq. (|23j) . as functions of a p , for a system with N = 768 amphiphiles. 
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Figure 9: Plot of the measured surface tension £ and of the estimate of the same quantity, as given 
by eq. l|23|) . as functions of a p , for a system with N = 1152 amphiphiles. 
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Figure 10: Plot of the effective fluctuation modes N' as a function of the vanishing tension projected 
area A* (a = 0) and as a function of the number of molecules N (upper x-axis). The line is a 
linear fit of the data, whose slope, as a function of A*, is 0.086 ± 0.003 ~ l/(47r), as predicted by 
eq. (JH. 
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